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Controlled aggregation of magnetic ions in a semiconductor. 
Experimental demonstration 

A. Bonanni,^'Ll A. Navarro- Quezada/ Tian Li,^ M. Wegscheider/ Z. Matej,^ 
V. Holy,2 R.T. Lechner,^ G. Bauer,^ M. Kiecana,^ M. Sawicki,^ and T. Dietl^'"^ 

^Institute for Semiconductor and Solid State Physics, 

Johannes Kepler University, A-4O4O, Linz, Austria 

^Department of Condensed Matter Physics, Charles University, Cz-121 16 Prague, Czech Republic 

^Institute of Physics, Polish Academy of Sciences, PL-02-668 Warszawa, Poland 

"^ Institute of Theoretical Physics, University of Warsaw, PL- 00- 681 Warszawa, Poland 

The control on the distribution of magnetic ions into a semiconducting host is crucial for the 
functionality of magnetically doped semiconductors. Through a structural analysis at the nanoscale, 
we give experimental evidence that the aggregation of Fe ions in (Ga,Fe)N and consequently the 
magnetic response of the material are affected by growth rate and co-doping with shallow impurities. 
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There is an increasing amount of evidence that owing 
to specific features of magnetic impurities in wide band- 
gap semiconductors and oxides, the epitaxial growth of 
these systems can result in the self-organized aggregation 
of magnetically robust nanocrystals embedded in the host 
paramagnetic matrix [l|, l2|, |3[ . With no doubt this finding 
holds enormous potential for the fabrication of a range of 
multifunctional nanosy stems relevant to spintronics, na- 
noelectronics, photonics, and plasmonics [2, 4]. However, 
it has also been realized [5] that enduring difficulties in 
the experimental resolution and identification of the em- 
bedded nanostructures hamper the progress in the visu- 
alization, understanding, and control of the mechanisms 
accounting for relevant and hitherto unexplored nano- 
assembly processes. 

In this Letter, by exploiting state-of-the-art nano- 
characterization tools, we show how growth conditions 
and co-doping with shallow donors or acceptors affect 
the distribution of the magnetic ions in (Ga,Fe)N de- 
posited by metalorganic vapor phase epitaxy (MOVPE), 
emphasizing conclusions that are pertinent to the whole 
class of wide band-gap diluted magnetic semiconductors 
(DMS) and diluted magnetic oxides (DM0), extensively 
studied over the last eight years [jj, ly, lZ|. In particu- 
lar, we identify different ways by which transition metal 
(TM) impurities can incorporate into the host semicon- 
ducting lattice and we link the structure investigated at 
the nano-scale to the macroscopic magnetic properties. 
Our findings show that Fermi-level tuning by co-doping 
with shallow impurities is instrumental in controlling the 
magnetic ions aggregation. This provides an experimen- 
tal support to recent theoretical suggestions [8|, |9| and 
it demonstrates that the Fermi-level engineering docu- 
mented so-far for iodine and nitrogen doped (Zn,Cr)Te [5] 
operates also in the case of III- V DMS. Moreover, we find 
that the aggregation of magnetic ions can be hindered by 
increasing the growth rate. This indicates that the nano- 
assembling process is controlled by a kinetic barrier to 



the surface diffusion of the magnetic ions, whose height 
depends on co-doping. 

All studied epilayers have been fabricated by MOVPE 
on c-plane sapphire substrates employing the growth pro- 
tocols and the in situ and ex situ characterization meth- 
ods we have reported previously for (Ga,Fe)N [10, [ll] 
and GaN:Mg [12], and that we apply now to produce 
(Ga,Fe)N and (Ga,Fe)N:Si,Mg layers. The total Fe 
concentration in the samples varies from 4x10^^ cm~^ 
to 3x10^^ cm~^ for the Fe-precursor (Cp2Fe) flow-rate 
ranging from 50 to 300 standard cubic centimeters per 
minute (seem) [10]. The Mg-concentration is found to be 
in the range 2-3x10^^ cm~^ [12] and the Si-content is 
estimated to be 1x10^^ cm~^. The growth-rate during 
the deposition is regulated by the Ga-precursor (TMGa) 
flow-rate and varies about linearly from 0.2 to 0.3 nm s~^ 
for 5 to 15 seem of TMGa flow, respectively. 

The difference between the magnetization values mea- 
sured by SQUID magnetometery up to 5 T at 1.8 K and 
5 K is employed to determine the concentration npg3+ of 
paramagnetic Fe^+ ions in the layers [ll|l, which agrees 
within a factor of two with the value of n-p^3+ obtained 
from electron spin resonance measurements [10]. Low- 
field and low-temperature hysteresis loops allow us to 
evaluate the concentration of Fe contributing to the ferro- 
magnetic signatures, nperro, which we determine assum- 
ing a magnetic moment of 2jj.b per Fe ion. 

In addition to high-resolution transmission elec- 
tron microscopy (HRTEM) and energy dispersive spec- 
troscopy (EDS) [10], we fiave performed and present fiere 
Fourier-filtering of the TEM images [13] for strain analy- 
sis [14], reconstructed by using the 002 spatial frequency. 
Powder x-ray diffraction (XRD) measurements [15], have 
been carried out at the beamline ID31 of the ESRF 
(Grenoble - France) using a photon energy of 15.5 keV. 
The x-ray data correspond to the diffracted intensities in 
reciprocal space along the sample surface normals, col- 
lected using a secondary crystal monochromator in front 
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FIG. 1: (color online) Synchrotron radiation powder XRD 
spectra vs. Fe content and growth rate for (Ga,Fe)N, as indi- 
cated. The two uppermost traces are generated by the same 
sample, but measured at different rotation angles (therefore 
different tilt). Inset: normalized number of Fe ions contribut- 
ing to the ferromagnetic signatures (dots) and near-band-edge 
PL intensity (triangles) vs. TMGa flow-rate, i.e. growth-rate. 



of the detector, thus the GaKc^ fluorescence is suppressed. 
We have acquired symmetric uj /2Q scans for all samples 
with scattering angles 26 up to 150 deg and found that 
the most intense peaks of the detected crystallographic 
phases appear in the 20-range up to 25 deg. As at the 
powder diffraction beamline the sample tilt can not be 
adjusted, neither the position nor the intensity of the 
substrate peaks may be compared for different scans (see 
the difference between the two upmost curves in Fig. 1). 
In contrast, the peaks for small nanocrystals are broad, 
and thus virtually insensitive to the sample tilt. 

It has already been established that once the solubil- 
ity limit of Fe in GaN at the given growth conditions 
is exceeded, the ferromagnetic response of (Ga,Fe)N in- 
creases with the Cp2Fe flow-rate [10|. In Fig. 1 we report 
the 00 /2Q powder diffraction synchrotron XRD scans for 
(Ga,Fe)N samples with different Fe-content above its sol- 
ubility limit into the GaN host (Cp2Fe 250 seem and 
300 seem, respectively). For the samples deposited at 
the low growth-rate (TMGa at 5 seem), in contrast to 
laboratory high-resolution XRD which does not evidence 
any phase separation in (Ga,Fe)N [10], the high intensity 
available at the synchrotron beamline allows to reveal the 
presence of new diffraction peaks identified as the (002) 
and (111) of the phase e-FesN (hexagonal siderazot struc- 
ture, space group No. 182 (P6322), lattice parameters 
a = 0.4698 nm and c = 0.4379 nm. Curie temperature 
Tc = 575 K and magnetic moment per Fe ion m = 2iib 
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This assignment is consistent with the crystallo- 



graphic characteristics of the Fe-rich secondary phases as 
put on view by TEM images, and summarized in Fig. 2. 
The height of the e-FeaN-related diffraction maxima is 
enhanced with increasing the flow-rate of Cp2Fe, how- 
ever no significant dependence of their full-width-at-half- 
maxima (FWHM) is observed upon varying the nominal 
content of magnetic ions. This implies that the mean 
size of the precipitates does not substantially vary with 
increasing nominal Fe content, whereas the density of the 
precipitates is enhanced, as confirmed by TEM [10]. We 
are, thus, led to the conclusion that the precipitation oc- 
curs by a nucleation mechanism in which only nanocrys- 
tals with a critical size can form. We use the FWHM 
of all diffraction maxima of e-FesN for an estimate of 
the mean size of the precipitates and by employing the 
Williamson-Hall plot method [l7] we obtain an average 
value for the nanocrystals diameter of 15 ± 5 nm. 

Our studies for various growth rates demostrate that 
the nanocrystal nucleation is limited by a kinetic bar- 
rier for the surface Fe diffusion. As shown in Fig. 1 the 
diffraction peaks originating from the Fe-rich nanocrys- 
tals appear only in the slow growth-rate regime (samples 
grown with TMGa at 5 seem) and are quenched at higher 
growth-rates, as evidenced by the lowest XRD trace in 
Fig. 1, obtained for the layer deposited at 10 seem TMGa. 
The inset of Fig. 1 gives the normalized number of Fe- 
lons contributing to the ferromagnetic signatures, nperro, 
as a function of the TMGa precursor flow-rate {i.e. the 
growth-rate) and its decrease with increasing flow-rate 
agrees with the XRD data, giving evidence of a reduced 
contribution of Fe-rich magnetic nanocrystals at higher 
growth-rates. This is further supported by a dramatic 
drop of the near-band-edge (3.49 eV at 10 K) photolu- 
minescence (PL) intensity, also reported in the inset to 
Fig. 1, witnessing, with increasing growth-rate, the en- 
hanced incorporation of the magnetic ions in unpaired 
positions, where they act as PL inhibitors [6]. 

At intermediate values of growth rate {e.g. TMGa 
at 10 seem) and a Fe-precursor flow-rate of 300 seem, 
where with synchrotron XRD the onset of second phases 
is not detected, ferromagnetic signatures are still clearly 
seen by SQUID up to a blocking temperature typically 
over 300 K. These puzzling observations are elucidated by 
the TEM images with mass contrast and strain mapping, 
shown in Figs. 2(g-l), that put into evidence the presence 
of spinodal decomposition into nano-scale regions with 
high Fe content embedded in the Fe-poor matrix with- 
out any crystallographic phase separation. This appears 
as a generic property of a number of DMS and DM0, 
in which no precipitates are detected, but ferromagnetic 
features persist up to high temperatures [7]. 

Remarkably, our TEM, XRD, and SQUID data reveal 
that the aggregation of Fe ions can be diminished or even 
prevented by co-doping with either Si donors or Mg ac- 
ceptors. Figure 2 presents TEM data for the two rel- 
evant initial regimes, namely (Ga,Fe)N with embedded 
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FIG. 2: (color online) Bright-field images (a, g), HRTEM with mass contrast (b, d, e, h, j, k) and Fourier filtered images with 
strain mapping (c, f, i, 1) of (Ga,Fe)N revealing the presence of FesN precipitates (a, b, c, e, f), spinodal decomposition (g-1), 
and the effect of co-doping by either Si (d, j) or Mg (e, f, k, 1) preventing the formation of the Fe-rich regions. 



Fe-rich nanocrystals [Figs. 2(a-f)] evidenced by Moire 
fringes contrast and (Ga,Fe)N showing spinodal decom- 
position [Figs. 2 (g-1)] generating mass contrast and lat- 
tice distortion, as proved by the filtered images [Figs. 2(i, 
1)]. From Figs. 2(d-f) and 2(j-l) the reduced aggregation 
of Fe-rich regions as a consequence of co-doping with Si 
and Mg, respectively, is evident. 

This effect is further corroborated by the XRD results 
given in Fig. 3. In both cases, the shallow impurities 
are found to hamper efficiently the precipitate aggrega- 
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FIG. 3: (color online) Effect of Si and Mg co-doping on the 
synchrotron radiation powder XRD spectra of (Ga,Fe)N; in- 
sets: effect of co-doping (a) with Si and (b) with Mg on the 
normalized number of Fe ions contributing to the ferromag- 
netic signatures in (Ga,Fe)N vs. TMGa flow-rate, i. e. growth- 
rate (a) and vs. Cp2Fe flow-rate (b). 



tion, so that the diffraction peaks corresponding to e- 
FesN are suppressed in the case of co-doped samples. 
We point out that the XRD curve from the Mg co-doped 
sample exhibits also a broad maximum identified as the 
(111) diffraction from pure 7-Fe (austenite). This finding 
is supported by our TEM observations, revealing in the 
considered sample the presence of a few-nm thick Fe in- 
clusions at the surface. The suppression of the ferromag- 
netic contribution in the co-doped layers is further vali- 
dated by the reduced number of average Fe ions adding 
to the ferromagnetic response in both Si- and Mg-co- 
doped (Ga,Fe)N layers. In inset a) to Fig. 3, the normal- 
ized nperro IS givcu for a constant Cp2Fe as a function of 
the TMGa flow-rate. As seen here, the quenching of the 
ferromagnetic contribution by Si co-doping is equally ob- 
served in samples deposited at a low growth-rate (5 seem 
TMGa, i. e. presenting Fe-rich nanocrystals when not co- 
doped) and at intermediate growth-rate (10 seem TMGa, 
i. e. showing - when not co-doped - spinodal decompo- 
sition). The change in the normalized nperro upon co- 
doping with acceptors is presented in inset b) as a func- 
tion of the nominal Fe content in the low magnetic ions 
doping regime (Cp2Fe at 50-150 seem). 

In order to explain these key findings, we note that ex- 
cept for Mn in II- VI compounds [5| , owing to the presence 
of the open d shells in the vicinity of the Fermi level, the 
nearest neighbor pair of TMs in semiconductors shows a 
large binding energy which promotes the magnetic ions 
aggregation [2|, |5|, [iS^ . However, if carriers introduced by 
co-doping can be trapped by these ions, the pair binding 
energy will be altered, usually reduced by the correspond- 
ing Coulomb repulsion [5|, [sl, l9|. While the presence of 
the mid-gap electron trap, z. e., the Fe+^/Fe+^ state, is 
well established in GaN 0, the level Fe+VFe+'^ is ex- 
pected to reside rather in the valence band [6']. However, 
it has been recently suggested that in GaN the potential 



introduced by the Fe+^ ion is strong enough to trap a 
hole in a Zhang- Rice type of state [19]. Therefore, since 
the concentrations of Si and Mg are comparable to that 
of nperro l^ our sauiplcs, the obstructive effect of Si and 
Mg co-doping on the nanocrystal formation and, thus, 
on the ferromagnetic response, is elucidated. Moreover, 
in order to test the scenario that the Fe aggregation oc- 
curs at the growth surface we deposited films in which 
Fe and Mg have been alternatively supplied in the ^-like 
fashion. No influence of co-doping on the ferromagnetism 
has been found, despite that Mg is known to diffuse in 
GaN [12]. Significantly, in view of our results, the previ- 
ously observed effect of co-doping on ferromagnetism in 
(Ga,Mn)N, and assigned to the dependence of the double 
exchange mechanisms of the spin-spin coupling on the po- 
sition of the Fermi level with respect to the center of the 
d band [20|, has to be reconsidered. Furthermore, since 
with the experimental resolution of 5x10^^ ijlb cm~^ we 
do not observe any spontaneous magnetization in GaN, 
GaN: Si and GaN:Mg down to 5 K, we conclude that the 
recently invoked vacancy-related ferromagnetism of GaN 
|2l| is not present in our samples. 

In summary, by combining HRTEM and synchrotron 
XRD with SQUID we have identified three distinct ways 
by which Fe incorporates into the GaN lattice: (i) sub- 
stitutional Fe^+ diluted ions accounting for the para- 
magnetic response [10|; (ii) Fe-rich (Ga,Fe)N wurtzite 
nanocrystals commensurate with and stabilized by the 
GaN host lattice and (iii) hexagonal e-FesN precipitates. 
The formation of nanocrystals containing a large density 
of the magnetic constituent elucidates the origin of the 
ferromagnetic features persisting up to above room tem- 
perature. Importantly, the co-doping with either Si or 
Mg hampers the nanocrystal assembling. This demon- 
strates that the charging of the magnetic ions, Fe^+ -^ 
Fe^+ and Fe^+ -^ Fe^+ + /i, respectively, inhibits the Fe 
aggregation and explains the sensitivity of the ferromag- 
netic response to co-doping with shallow donors or accep- 
tors. Furthermore, the influence of the growth rate on the 
nanocrystal formation indicates that the Fe aggregation 
occurs at the growth surface. Since, quite generally, the 
binding energy of TM pairs depends on the valency of the 
open (i-shells, there is a ground to suppose that the Fermi 
level engineering evoked here for (Ga,Fe)N:Si,Mg can 
serve to control the magnetic ion aggregation in a num- 
ber of semiconductors and oxides, providing a way to the 
self-organized fabrication of multi-component systems 
with tailored magnetic, magneto-optical, and magneto- 
transport properties at the nanoscale. 
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